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At the Fermilab Tevatron energies, (-^=1800 GeV and ^=630 GeV), jet production is the dominant process. 
During the period 1992-1996, the D0 and CDF experiments accumulated almost 100 pb _1 of data and performed 
the most accurate jet production measurements up to this date. These measurements and the NLO-QCD theo- 
retical predictions calculated during the last decade, have improved our understanding of QCD, our knowledge of 
the proton structure, and pushed the limit to the scale associated with quark compositeness to 2.4-2.7 TeV. In this 
paper, we present the most recent published and preliminary measurements on jet production and fragmentation 
by the D0 and CDF collaborations. 



1. Introduction 

Quantum chromodynamics (QCD) describes 
the inelastic scattering between a proton and an 
antiproton as a hard collision between their con- 
stituents or partons: quarks or gluons. After 
the collision, the outgoing partons hadronize into 
streams of particles called jets. At the Tevatron 
energies, (^=1800 GeV and yi=630 GeV), jet 
production is the dominant process. During the 
period 1992-1996, the D0 and CDF experiments 
accumulated almost 100 pb _1 of data and per- 
formed the most accurate jet production measure- 
ments up to this date. Among the results pub- 
lished in that period and subsequent years, we 
can cite inclusive jet cross sections, dijet angu- 
lar distributions, and dijet mass cross sections [ 
0, |, §, |§ § §, 0. At the same time, predic- 
tions for jet production rates have improved in the 
early nineties with next-to-leading order (NLO) 
perturbative QCD calculations [ [| and more ac- 
curate parton distribution functions (pdf) [ |To| . 
The high center-of-mass energy at the Tevatron 
and the unprecedented accuracy of the measure- 
ments, together with the NLO-QCD theoretical 
predictions derived during the last decade, have 
improved our understanding of QCD, our knowl- 
edge of the proton structure, and pushed the limit 
to the scale associated with quark compositeness. 

In this paper, we include a summary of some of 
the most significant jet results published by CDF 
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and D0 as well as their most recent preliminary 
measurements. Jet cross sections in forward pseu- 
dorapidity regions, cross sections of dijets sepa- 
rated by large pseudorapidity intervals, and sub- 
jet and particle multiplicity measurements pro- 
vide information on parton distribution functions, 
probe BFKL dynamics, explore the jet structure, 
and study the hadronization process. In the first 
sections, we describe how jets are selected, recon- 
structed and calibrated at D0 and CDF. They 
are followed by sections on each measurement and 
the conclusion. 

2. Jet Reconstruction and Data Selection 

For most of the analyses presented here jets are 
reconstructed using an iterative fixed cone algo- 
rithm with a cone radius of 1Z=0.7 in rj-cj) space [ 
[ill , (pseudorapidity is defined as 77 = — ln[tan|]). 
This algorithm is applied to calorimeter towers 
without making use of tracking information, ex- 
cept for the determination of the interaction ver- 
tex. The D0 subjet multiplicity measurement 
uses a Kt algorithm [ |24|, Q on calorime- 
ter towers, with a resolution parameter T>=1 (see 
Ref. [ 0.) The CDF particle multiplicity re- 
sults use fixed cone algorithms with different cone 
sizes, based on particle information at the track- 
ing level. 

The offline data selection procedure eliminates 
background caused by electrons, photons, noise, 
or cosmic rays. In the case of D0 it follows the 



methods described in Refs. [ [H| |TJ]. 

3. Energy Corrections 

The jet energy scale correction, described by 
D0 in Ref. [ [l5|], removes instrumentation effects 
associated with calorimeter response, showering, 
and noise, as well as the contribution from spec- 
tator partons (underlying event). 

The D0 energy scale correction corrects the jet 
Et from their reconstructed value to their "true" 
Et on average (energy of a jet defined from fi- 
nal state hadrons). An unsmearing correction is 
applied later to remove the effect of a finite Et 
resolution [ 13 1. CDF corrects both for scale and 



resolutions using a Monte Carlo simulation tuned 
to represent the data. 

4. The Inclusive Jet Cross Section at 
= 1800 GeV 

The inclusive jet cross section is measured by 
both the D0 (in \rj\ < 0.5 and 0.1 < \rj\ < 0.7) 
and the CDF experiments (in 0.1 < \r]\ < 0.7). It 
is defined as: 

dErdr] deiSET^-r] 

where Ni is the number of accepted jets in Et bin 
i of width AEt, C-i is the integrated luminosity, 
ti is the efficiency of the trigger, vertex selection, 
and the jet quality cuts, and Ary is the width of 
the pseudorapidity bin. 

Figure [l] shows the various uncertainties for 
the D0 (\r)\ < 0.5) cross section. The second 
outermost curve shows the error on the energy 
scale which varies from 8% at low Et to 30% at 
450 GeV and dominates the total error. Most of 
the systematic uncertainties of the inclusive jet 
cross section are highly correlated as a function 
of Et- 

Figure || show the fractional difference between 
the data, D, and a jetrad theoretical prediction, 
T, normalized by the prediction, ((D — T)/T), for 
\r)\ < 0.5. The jetrad prediction was generated 
with [i = 0.5ETmaxi Ttscp — 1-3 and several dif- 
ferent choices of pdf. The error bars represent 
statistical errors only. The outer bands represent 
the total cross section error excluding the 5.0% 
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Figure 1. Contributions to the D0 jet (\r]\ < 0.5) 
cross section uncertainty plotted by component. 



luminosity uncertainty. Given the experimental 
and theoretical uncertainties, the predictions are 
in agreement with the data; in particular, the 
data above Et = 350 GeV show no indication 
of an excess relative to QCD. 

The data and theory can be compared quanti- 
tatively with a x 2 test incorporating the uncer- 
tainty covariance matrix. The \ 2 is given by: 



X 



(2) 



1,3 



where Si is the difference between the data and 
theory for a given Et bin, and is element i, j 
of the covariance matrix: 



ViA = p,, ■ Aov. • Act,- 



(3) 



where Act is the sum of the systematic error and 
the statistical error added in quadrature if i = j 
and the systematic error if i ^ j, and pij is the 
correlation between the systematic uncertainties 
of Et bins. 

All but one of the jetrad predictions ade- 
quately describe the \rj\ < 0.5 and 0.1 < |??| < 0.7 
(shown in Figure ||) cross sections (probabilities 
for x 2 to exceed the calculated value are between 
10% and 86%). The prediction using CTEQ4HJ 
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Figure 2. The difference between D0 data and JE- 
trad QCD predictions normalized to predictions 
for \rj\ < 0.5. The shaded region represents the 
±lcr systematic uncertainties about the prediction. 

and fi = 0.5E produces the highest probabil- 
ity for both measurements. The prediction with 
the MRSTGD pdf has a probability of agreement 
with the data of 0.3%, and is incompatible with 
the data. 

The top panel in Fig. | shows (D - T)/T for 
the D0 data in the 0.1 < \r/\ < 0.7 region rela- 
tive to a jetrad calculation using the CTEQ4HJ 
pdf, /i = 0.5ETmax, and 72 SO p = 2.072.. Also in- 
cluded is the published CDF measurement from 
the 1992-93 Tevatron running period [ |^| relative 
to the same JETRAD prediction. The CDF mea- 
surement shows an excess with respect to the the- 
ory at high Et, which can be accommodated by 
adjusting the gluon pdf (CTEQ4HJ set). If we 
include the systematic uncertainties of the two 
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Figure 3. Top: Normalized comparisons of the D0 
data and of the CDF data to a jetrad predic- 
tion (CTEQ4HJ and /i = 0.5E Tmax ). Middle: Dif- 
ference between the data and smoothed results of 
CDF normalized to the latter. The shaded region 
represents the ±1ct systematic uncertainties about 
the D0 data. The dashed curves show the ±1<7 
systematic uncertainties about the smoothed CDF 
data. Bottom: A comparison of the systematic un- 
certainties of both experiments. 

experiments (CDF's uncertainties in Ref. [ ^ ) in 
a covariance matrix, the x 2 is 30.8 for 24 degrees 
of freedom (probability of 16%), representing ac- 
ceptable agreement between D0 and CDF. 

5. Inclusive Jet Cross Section at \/s— 
630 GeV 

Figure || shows the fractional difference be- 
tween the D0 data and several jetrad predic- 
tions given different choices of renormalization 
scale and pdf. These NLO QCD predictions are 
in reasonable agreement with the data. The data 
and predictions are compared quantitatively with 
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Figure 4. (D-T)/T for the D0 630 GeV jet cross 
section. (\r)j e t\ < 0.5). The solid stars represent the 
data compared to the calculation for /i = O-BExmax 
and the pdfs CTEQ4M , CTEQ4HJ , MRST , 
MRSTGU and MRSTGD . The shaded region rep- 
resents the ±lcr systematic uncertainty, 
a x 2 test- All but two of the JETRAD predictions 
adequately describe the cross section at *Js = 
630 GeV (the probabilities for x 2 to exceed the 
calculated values are between 10% and 74%). The 
prediction using MRSTGU and [i = 0.5ETmax 
produces the highest probability. The predic- 
tion with MRSTGD pdf and (jt = 0.5E Tmax , and 
CTEQ3M pdf and fi — 2ETmax are ruled out 
by the D0 measurement (agreement probability 
<0.4%). 

6. The Ratio of Jet Cross Sections 

The dimensionless inclusive jet cross section is 
given by 



& CM 



d 2 a 
2n dE^drj ' 



(4) 



where d 2 a/dExdr] is given by Eq. |l|, and x is the 
center-of-mass energy. A naive parton model with 
no Q 2 dependence of the pdfs, and therefore no 
running of a s , predicts this cross section ratio 
to be unity, that is independent of the center-of- 
mass energy of the pp system. The observable is 
nearly insensitive to the choice of parton distribu- 
tion functions. It is, therefore, a more stringent 
test of QCD matrix elements. 

D0 has measured the ratio of inclusive jet 
cross sections in the central pseudorapidity bin 
<0.5) [||. This quantity is calculated in bins 
of identical xt- 



R(x t ) 



C630 (x T ) 
C1800 ( x t) 



(5) 



The inclusive jet cross section errors are highly 
correlated as a function of Et and center-of-mass 
energy and will cancel in the ratio. The energy 
scale uncertainty dominates the total error in the 
ratio too. 

Figure || shows the ratios of cross sections 
with jetrad predictions using different pdfs. 
The measured ratios lie approximately 10% be- 
low the theoretical predictions, which have an 
uncertainty of about 10%. The x 2 values lie 
in the range 15.1-24 for 20 degrees of freedom 
(corresponding to probabilities in the range 28% 
to 77%). The best agreement occurs for ex- 
treme choices of renormalization scales (/x = 

0.25, 2ETmax)- 

In general, the NLO-QCD predictions yield sat- 
isfactory agreement with the D0 data for stan- 
dard choices of renormalization scale or pdfs. In 
terms of the normalization, however, the absolute 
values of the standard predictions lie consistently 
and significantly higher than the data. 

CDF has also measured the ratio of central jet 
cross sections. Figure ^ shows a preliminary re- 
sult by CDF compared with NLO-QCD predic- 
tions (and the D0 measurement). Although the 
data and the theory agree in shape at high xt, 
there is a significant deviation at low values. 
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Figure 5. The D0 ratio of dimensionless 
cross sections compared with jetrad predictions 
= 0.5E Tmax and the CTEQ3M , CTEQ4M , 
CTEQ4HJ , or MRST pdfs). The shaded band 
represents the ±lcr systematic uncertainty. 
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7. The Forward Jet Cross Sections at 
y/s = 1800 GeV 

D0 has performed preliminary measurements 
of forward jet cross sections up to pseudorapidi- 
ties of \rj\ = 3. These measurements allow to 
reach regions in (x,Q 2 ) space previously unex- 
plored. Figure [?] shows a comparison of the mea- 
sured cross sections in five different r\ bins with 
the jetrad prediction (CTEQ4HJ or MRST, and 
/i = ETmax/2)- Within the experimental and the- 
oretical uncertainties, the measurement and the 
calculation are in good agreement, x 2 studies are 
underway for a quantitative statement. 
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Figure 7. Pseudorapidity dependence of the in- 
clusive jet cross section {\rj\ <3), compared with 
NLO QCD (CTEQ4HJ: full circles, MRTS: open 
circles). The bands represent the total system- 
atic uncertainties in the experiment. 



Figure 6. The D0 and CDF ratio of dimensionless 
cross sections compared with each other and the 
jetrad predictions. 



8. The Ratio of Dijet Mass Spectrums at 
y/s = 1800 GeV 

The dijet mass spectrum is calculated using the 
relation: 



d 3 a 



dMjjdrjidr]2 C i e l AMjjArnAr]2 



(G) 



where 2Vj is the number of events in mass bin 
i; Ci is the integrated luminosity; e{ is the effi- 
ciency of the trigger, vertex selection, and the jet 
quality cuts; AMjj is the width of the mass bin; 
and Ar/i t 2 are the widths of the pseudorapidity 
bin At D0 the cross section is measured for the 
pseudorapidity bin \rjj et \ < 1.0 . The systematic 
errors are dominated by the uncertainties due to 
the jet energy scale, which are 7% (30%) for the 
209 (873) GeV mass bins. 

The dijet mass cross section measurement was 
then repeated for \r]j et \ < 0.5, and 0.5 < \rfj e t\ < 
1.0 and their ratio was determined. A large frac- 
tion of the total error cancels, as well as the un- 
certainty in the theoretical prediction of the ratio 
which is less than 3% due to the choice of pdf, and 
6% from the choice of renormalization and factor- 
ization scale (excluding /i = 0.25ETmax)- By tak- 
ing the ratio n (\r]jet \ < 0.5) /k (0.5 < \r)j e t\ < 1-0) 
the systematic uncertainties decrease to less than 
10%. 

Given the experimental and theoretical uncer- 
tainties, the prediction can be regarded as in good 
agreement with the data (see Fig. |J). The data 
are also in agreement, within the uncertainties, 
with the cross section measured by CDF [ ^| . 

All choices of pdfs and renormalization scales 
are in good agreement with the data (x 2 ) test, 
except for /j = 0.25ETmax which is excluded by 
the data. 

The ratio of the mass spectra is used to place 
limits on quark compositcncss. The pythia event 
generator is used to simulate the effect of compos- 
iteness by taking the ratio of these LO predictions 
with compositeness, to the LO with no compos- 
iteness, and scaling with this factor the jetrad 
NLO prediction (shown in Fig. ||). 

The D0 data shows no evidence of compos- 
iteness. The dijet mass spectrum rules out quark 
compositeness models at the 95% confidence level 



06 



in 

no 0.75 



u 



D0 Data 
A + =1.5TeV 

LL 

A + =2.0TeV 

LL 

■ A + =2.5TeV 

LL 

• A + =3.0TeV 

LL 




M (GeV/c ) 

Figure 8. The ratio of dijet mass cross sections for 
\r)j et \ < 0.5 and 0.5 < \r)j et \ < 1.0 for data (solid 
circles) and theoretical predictions for composite- 
ness models with various values of Aj L . The error 
bars show the statistical uncertainties. The shaded 
region represents the ±1ct systematic uncertainty. 



where Aj L is below 2.7 TeV and A LL is below 
2.4 TeV. 



9. The Triple Differential Jet Cross Sec- 
tions at yfs — 1800 GeV 

Both CDF and D0 have performed preliminary 
measurements of triple differential jet cross sec- 
tions. This observable is defined only in terms 
of the two leading jets of the event, to provide 
information on pdfs. In particular, the Tevatron 
probes high Q 2 and x values previously unreach- 
able. 

CDF defines the triple differential jet cross sec- 
tion as dET f d ^ ldr)2 ' where E T \ is the transverse 
energy of a central jet in the event, and 771 , rfl are 



the pseudorapidities of the central and forward 
jets, respectively. This quantity is measured as 
a function of the transverse energy of the central 
jet for different 772 bins up to |t?|=3. 



D0 measures 



versus Et , where Et is 



dETdrj\drj2 

the transverse energy of the central or the forward 
jet (the event enters twice in the measurement). 
This quantity is measured for different 772 bins up 
to 1 77 1 =2. The main difference between the CDF 
and D0 observables is that D0 measures the Et 
of both jets and CDF measures only the Et of 
the central jet. 



(Data-Theory )/Theory, Theory = CTEQ4M 
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Figure 9. D0 triple differential jet cross sections in 
different \r]2\ < 2 intervals for the two leading jets 
(central and forward) in the same pseudorapidity 
side (sign of 771 same as sign of 772). 
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Figure 10. D0 triple differential jet cross sections 
in different \r]2\ < 2 intervals for the two leading jets 
(central and forward) in opposite pseudorapidity 
sides (sign of 771 opposite of sign of 772). 



Figures |p- [h] show the D0 measurement in dif- 
ferent 772 bins for the two leading jets (central and 
forward) in the same pseudorapidity side (sign of 
771 same as sign of 772), and opposite sides. Each 
configuration is adequate for learning about dif- 
ferent pdfs in different regions of (x,Q 2 ) space. 
For example, two forward jets in the same side are 
associated with one incoming parton with high x 
and the other with low x. Figure [ll] shows the 
CDF result for different 772 bins, independently of 
the relative sign of 771 and 772. 

Qualitatively, there is good agreement between 
data and theory. Both experiments are currently 
working on quantitative studies. 

10. Dijet Cross Sections at Large 7; Inter- 
vals 

At high center-of-mass energies, y/s, and for 
momentum transfers, Q, fixed and <C \fs, the 
radiative corrections to the parton-parton scat- 
tering contain large logarithms ln(s/Q 2 ), which 
need to be summed to all orders in a s . This sum- 
mation is accomplished by the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) equation [ |l6||. 

Inclusive dijet production provides an ideal 
possible signature of BFKL dynamics. For large 



values of the jet longitudinal momentum fraction, 
Xj, the large logarithms ln(s/Q 2 ) result in large 
ln(s/Q 2 ) (where s/§ is the partonic center-of-mass 
energy) which factorize in the partonic dijet cross 
section, a. The \n(s/Q 2 ) terms are of the order of 
the pseudorapidity interval, A77, between the two 
jets [ |l7| (77 = — ln(tan(0/2)), where 9 is the polar 
angle of the jet relative to the proton beam). 

D0 performed a measurement of the dijet cross 
section at two different center-of-mass energies, 
y/sj, = 1800 GeV and = 630 GeV, using 

the D0 detector at the Fermilab Tevatron. The 
kinematics of the event is reconstructed using the 
most forward/backward jets, and the cross sec- 
tion is measured as a function of x\, xi and Q 2 
at each center-of-mass energy. The ratio of the 
cross sections is then determined at the same val- 
ues of xi, X2 and Q 2 between the two energies. 
This eliminates the dependence of the cross sec- 
tion on the pdfs and reduces the ratio to that of 
the partonic cross sections. It can be shown [ [?]] 
that the latter is a function only of the pseudora- 
pidity separations between the two jets (Ati's): 

&(Ar)A) e ( a BFKL-l)(At)A-A))B) 

R = "TTT r = ; = • (7) 

<7( A? 7s) y/ArjA/ArjB 
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Figure 11. CDF triple differential jet cross sections 
for different 772 bins. 

In other words, variation of y/s, while keeping x\, 
X2 and Q 2 fixed, is equivalent to variation of A77, 
which directly probes the BFKL dynamics. 

Several theoretical predictions can be com- 
pared to the D0 measurement. Leading Order 
QCD predicts the ratio of the cross sections to fall 
asymptotically toward unity. The herwig [ |T^| 
Monte Carlo provides a more realistic prediction. 
It calculates the exact 2-^2 subprocess includ- 
ing initial and final state radiation and angular 
ordering of the emitted partons. The LLA BFKL 
intercept for a s (20GeV) = 0.17 [ [IJ) is equal to 
1.45. The Next-to-Leading Logarithmic [ |2Q | are 
not as yet available. 

The ratio of cross sections is shown in Fig. |l2| 
as a function of the mean pseudorapidity interval 
at 630 GeV. It is evident that the growth of the 
dijet cross section with Aij is stronger in the data 
than in any theoretical model which was consid- 
ered. Namely, the measured ratio is higher by 
4 standard deviations than the LO prediction, 3 
deviations than the herwig prediction, and 2.3 
deviations than the LLA BFKL one. 



11. Sufojets»M Multiplicities 
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A jet is typically associated with the energy 
and momentum of each final state parton. Ex- 
perimentally, however, it is a cluster of energy in 
the calorimeter. QCD predicts that gluons radi- 
ate more than quarks. Asymptotically, the ratio 
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Figure 12. The D0 ratio of the dijet cross sec- 
tions at both center-of-mass energies for A77 > 1 
and A77 > 2. The inner error bars on the data 
points represent statistical uncertainties; the outer 
bars represent statistical and uncorrelated system- 
atic uncertainties added in quadrature. The error 
bars on the herwig predictions represent statisti- 
cal uncertainties. 

of objects within gluon jets to quark jets is ex- 
pected to be in the ratio of their color charges 
C A /C F = 9/A[^. 

D0 performed a preliminary measurement of 
subjet multiplicities in quark and gluon jets, as 
well as the ratio of the means of these two quan- 
tities. For this analysis, jets are reconstructed 
using the Kj- algorithm [ |2^, |2^, |24| with a reso- 
lution parameter T>=1 (see Rcf. [ |ll[). 

M is the subjet multiplicity in a mixed sample 
of quark and gluon jets. It may be written as a 
linear combination of subjet multiplicity in gluon 
and quark jets: 



M = fM g + (1 - f)M q 



(8) 



The coefficients are the fractions of gluon and 
quark jets in the sample, / and (1 — /), respec- 
tively. Consider Eq. (||) for two samples of jets 
in the same kinematic range, one at ^Js — 1800 
(gluon dominated) and the other at 630 GeV 
(quark dominated) , assuming M g and M q are in- 
dependent of y/s. The solutions are 



M q = 



Jl800jyj630 _ /630j^l800 
J1800 „ J630 



(9) 
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(10) 



where M 1800 and M 630 are the experimen- 
tal measurements in the mixed jet samples at 
y/s = 1800 and 630 GeV, and / 1800 and J 630 are 
the gluon jet fractions in the two samples. The 
method relies on knowledge of the two gluon jet 
fractions. 

Figure [l|| shows that the subjet multiplicity is 
clearly larger for gluon jets compared to quark 
jets. The gluon jet fractions are the largest 
source of systematic error. The measured ratio 
and its total uncertainty are: R = ^/ s j-i = 

1.91 ± 0.04(stat)+°;Jg(sys). The ratio is well de- 
scribed by the HERWIG parton shower Monte 
Carlo, and is only slightly smaller than the naive 
QCD prediction 9/4. 
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Figure 13. Corrected subjet multiplicity in quark 
and gluon jets, extracted from D0 data. 

12. Particle Multiplicities 

Perturbative QCD calculations, carried out in 
the framework of the Modified Leading Log Ap- 
proximation [ [2^] (MLLA), complemented with 
the Local Parton-Hadron Duality Hypothesis [ ^6) 
(LPHD), predict the shape of the momentum dis- 
tribution, as well as the total inclusive multiplic- 
ity, of particles in jets. The MLLA is an asymp- 
totic calculation, which proves to be infrared sta- 



ble, in the sense that the model cutoff parameter 
Qeff can be safely pushed down to Aqcd- LPHD 
is responsible for the hadronization stage and im- 
plies that hadronization is local and happens at 
the end of the parton shower development. In its 
simplest interpretation, the model has one param- 
eter Klphd, the rate of parton-to-hadron conver- 
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(11) 



In MLLA, momentum distributions and mul- 
tiplicities in quark and gluon jets in a restricted 
cone of size 9 around the jet axis are functions of 
Ejetd/Qef f [ and differ by a factor r: 

N<i-Mf£\ = Ijvs-ie*(£) 
r 



£ = log - , X = Ptrack/E je t 

X 



(12) 



Jets at the Tevatron are a mixture of quark and 
gluon jets. Therefore, 



n: 



charged 



hadrons 



(0 



IC charged/ . \ \ pnM LLA *j-q—jet fc\ _ 

^ lphd ytg y 1 g>~> iy part — 

KN^fiZ) (13) 

where e g is the fraction of gluon jets in the 
events, the factor of 1/r reflects the difference 
between gluon and quark jets, and, finally, the 
factor F nMLLA accounts for the next-to MLLA 
corrections to the gluon spectrum. Theoretical 
calculations [ Gq] predict somewhat different val- 
ues of F nML "\ but all agree that F nMLLA has 
almost no dependence on the jet energy in the re- 
gion relevant to this analysis. The average of the 
results above was chosen and the difference be- 
tween predictions was used as a theoretical error: 
F nMLLA =1.3±0.2. The same papers predict the 
value of r to be between 1.5 and 1.8. 

12.1. The Dijet Data Analysis 

CDF data collected during the 1993-1995 run- 
ning period was used for this analysis. Events 
with two jets well balanced in transverse energy 
were selected. Both jets were required to be in the 
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Figure 14. Fit of the parameter K for Klphd and 
r. Cone 0.47. First error is combined statistical 
and systematic errors, the second one - theoretical 
error coming from p nMLLA , 



central region. Tracks were counted in restricted 
cones of sizes 0.28,0.36 and 0.47 around the jet 
axis. 

Analysis of the fitted parameter K allows an 
extraction of both Klphd and r. According to 
Eq. [l3| the dependence is linear. Figure [l^, shows 
9 values of K (corresponding to 9 dijet masses for 
the largest cone-size 0.47) vs the gluon jet fraction 
(extracted using Herwig 5.6) in the events from 
respective dijet mass bins, as well as the results of 
the fit for Klphd and r. The same parameters 
can be extracted from the inclusive multiplicity 
using an integrated version of Eq. [13|. In this 
case, the extracted parameters will only rely on 
the total multiplicity and not on the exact shape 
of the distribution. Figure [l5| shows the fit of 
data with MLLA predictions as well as the fitted 
parameters Klphd and r. It is remarkable that 
the two results are in such a good agreement. 

12.2. Model-Independent Measurement 

The multiplicity in dijet and 7-jet events is 
compared (data selection was similar) to extract 
model-independent measurement of r. These 
samples have very different fraction of gluon jets 
for the jet energies 40-60 GeV (roughly 60% for 
dijets and 12% for 7-jet, according to Herwig 5.6). 
The multiplicities measured for each of the sam- 
ples and a knowledge of the gluon jet fractions 
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Figure 15. Charged particle multiplicity (per jet) 
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Figure 16. Ratio of charged multiplicities in gluon 
and quark jets based on comparison of the dijet 
and 7- jet events. 

allowed to extract r. Figure [l6] shows the mea- 
sured r as a function of the jet energy. The result 
for r is 1.75±0.11±0.15 in perfect agreement with 
MLLA result. 

13. Conclusions 

The 1992-1996 collider run at Fermilab repre- 
sented a major step in the testing of QCD. The 
D0 and CDF experiments measured jet cross sec- 
tions with unprecedented accuracy, extended the 
energy reach to ^450 GeV, and set a new limit 
of 2.4-2.7 TeV for the quark compositeness scale. 
In general, and within experimental and theoreti- 
cal uncertainties, QCD is in good agreement with 
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the data. Measurements on jet structure and 
fragmentation were also performed and yielded 
agreement with QCD; they support the perturba- 
tive nature of jet fragmentation. The upcoming 
run at the Tevatron, scheduled to start in March 
2001, will extend the energy frontier even further 
(y/s=2 TeV) and collect at least 20 times more 
data, allowing precision measurements of QCD 
in kinematic regions previously unexplored. 
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